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method using 2 telestope, high r2solution television camera, 
video tape recorder, d@geseal ~storage oscilloscopes, and 
calculator system has been devisad, tested and utilized. A 
jaser beam signal is recorded on video tape for further 
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Weerd a Tektronix 453 digital storage oscilloscope. The 
G@ig2tized signal is sant ¢5 2 Hewlett-Packard 9825 computing 


System for Fourier tcansform analysis and determination of 


@emeectractive index structure constant, rar Safeco! 
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ne BACKGROUND 

The increasing use of lasers and iaser technology for 
military applications has brought about 32 need for analysis 
Semeene laser beam in itS 2nvirorment, Pine (rtd bene 
atmosphere. A projerzt at Naval Postaqraduate Schooi dealing 
fee? S )§«6subject is continuinzs ard is the main topic of 


mies thesis. 
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The patterns produced by lasers on target 
roblems that includ2 broadening, beam wandér end intensitt 


Megetiaetions brought about by turbulence in the atmosphere, 


have been weli deternined [ Ref. 1]. DAG eccENSsOneene —=ac7d 
Moaet, {| Ref. 2], Cea Se ee ee ace ve index Pstrus ure 
Sonstant, has been determinei +o adequately express che 


Mem@meat.On Transfer Function (MIF), cr the Mutual Coherence 
meume@eson (MCF) for che atmospheres. 

A system has been developed that provides a measurement 
Or Cy? for atmospheric turbulznce along the optical path 


through which a laser is propajated. Tha system emplovs a 





Mmmeecon and t<elescops aS the datsctor and 2 distant laser as 
the source. This system duplicates the slit scanning system 
Peesertiy in use at Naval Postgraduete School [Ref. 2]. 
Tests using the vidicon equipment have been previously 
completed with values for C,? 2n the order of 10-15 n-2/3 
being obtained [Ref. 3 }. Measurements have been mad2 using 
the same experimental set-up, 2xcept that the Tektronix 468 
megetal Storage Jstilloscop>= are used for digitization 
¢nstead of the Quant2x DS-30 Dizital Video Processor. These 
feeeeieements have denonstrated values of 2, of comparable 


accuracy. 


B- PROBLEM 

Pemasing ~he digitel stoeraj: capability of the Tektronix 
468, the previously recorded laser sign2i is used 2S an 
input and is evaluat2ei by modifying the program developed by 
Crager [{Ref. 4]. A brieI overview of the precedure iis 
described below. Detailed explantions of the experiment and 
Sempwe St program ar2 contaionsi in Chapters [II and IV, 
respectively. 

The approach taken is basically the same as that of 
Crager, but because of the unavailability of a disk ROM it 


is necessary to store both raw and processed data on 
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Bime@me- iC  <2pe. Rieder ees (255 Loeided “and frecer iss - as 
necessary by the HP 3325. Ths basic asSunption vosrmitting 
M@eeysSsis is that th2 horizontal TV scan line through the 
Myset SpOt iS considered to accurately mirfcor 2 point spread 
function of the inage. Final analysis using the 468 
oscilloscope has shown this t> ba @ valid hypothesis. The 
Seeger ized data from the 468 oscilloscope agree with the 
previously measured lata which 152 many pixels of ~elevision 
ieee G2 g2tized by the DS-30 t95 express the point spread 
met iON. 

The sequence of analysis is that a video tape recording 
is made of the TV image of the propadcated laser bean. fe 
Semyeeut Of the vids recorder is sent <2 the Textronix 468, 
where the deérived TV scan lin2 is digitized, averaged, and 
Swered. The HP 9825 records th2 ligitized data and produces 
@a iine spreai function (LSF), DY -Gheegrating the ooint 


spread function, 


LSF(x) = | PSF(r) dy (e410) 


where 
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Smo computes ths Fourisr transform of ths LSF. The 
Meeeesaction limited Fouri2r cransform of the optics is now 
Sememeed if the MNIF sf the optics has not been previously 
measured experimentally. Next, the program finds the MIF of 


+he atmosphere by dividing the Fourier transform of the LS? 
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Seed, the program computes a single value for C,@. 


The program now predicts the size of 3 laser spot on 2 
Meeee Using the calsulated values of C,*. This value of C, 


is used t> calculate an MPF of th2 atmosphere which is chen 


Meee pitied by the Pourier transform of 2 source, 299  =ne 
Fourier transforn Dis} tn er el oe cS 5 Wiel = PlLOgds2m )) sneer 


Gameutates the invers2 Fourier transform of the products of 
the above and uses the Abel transform to give the angular 
Peemt spread intensity distriogution. Role ena sda aye oe 
fraction of energy as a functioa of the total energy within 


@ given radius R is calculated. 
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The thecry of laser beam propagation through 2 turbulent 
medium has been explained by Crittenden, aad others, and iis 
re-emphasized here for continuity purposes [Ref. 2]. S nee 
the effects of turbulence on laser weapons is of major 


concern, measurement and prediction capability for these 


VASemese ais 8s Cheers 7 ct aleyan 


O 
os 


ere Ccts 

Bae tO the existence of Tetiprocity, these ideas apply 
Memetth=r laser designators of imagers (Ref. 5]. ieee aoe 
stucy of the theoretical model by Fried, the effects of 
memeepheric turbuleizse ire ilivestigated [ Ref. 6}. This 
model uses euew sie Se 1a JONG Sen ODelCalm tlransce= *lune=i6n 
(OTF) when considering atmospnaearic turbulence. LDS ong 
meameOorTrr results from simply taxing an image of sufficiently 
long term which sees effectively all possibie turbulence 
Someaguratiocns. 

meen Observing ta= effect sf diffraction as the result 
Of the finite aperture of a point source, it can be seen 
that the resulting WMACee ie aoc Tune erm. COBS20eGCrsng he 
point spread function to have the same shape in the image 
plane regardless of its positisd2, the image function is «he 


result of eae sCOMVOLUELON of the source function and the 
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@@eecs diffraction function. [he convolution thneoren, as 


Gescribed in Fourier transfora theory, yields 


i (v (xX) ,v(Y)) =H (V(x) ov (Cy) } FO (V(X) ,V (y)) FM (V(X) Y (y)) (2 1) 


where 


POUL LOSER tEADororN Of -hev2nage function 


1 (v(x) -V(Y) ) 


eee) ,>V(Y)) = Fouriec transforn of the object function 


M(v(x),(v(y)) = @ Yes iter sFuUuees on OL 
: = 


v = svatial frequency 


Puepeaht SOUrcS such 3s 2 laser carn b2 analyzed in twe 
Peens7Ons using a point spreai function. Parough- Fourier 
Seenstorm theory, Phe. Stade Pe=ngP Spread. funceccion ss 
[eens rtOrMed into a two dimensional OTF of the optical 
Signal. This problem may be simplified by scanning the 
emege point spread function usiag a vidicon or slit-scanning 
system. Pniedoplvensecne Convi lution theogten, “the Fouztier 
meets tolm of the image point spread function is multiplied 
Meemeche Fourjer transform of the optical system resulting 


in the Fourier transform of th? overall system. The Abel 
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dimensional image PSF [Ref. 7]. 
As demonstrated by Crittend2n, and others, a numerical 
value for Cp? may b2 obtained by curve fitting using «he 


fol Ltowing: 

M = exp(-21.49*C ,2*Z *#£5 /3*\-1/3) (2.2) 
where 
M = MTF of the atmosohere 


Z= Range in meters 


ri 
if 


Py = Angular spatial frequency in cycles/radian 


ry 
ii 


rocal length of the optical system 

v = linear spatial fr2quency in cycies/meter 

A = wavelength in net2rs 

Cyne is obtained by a linear regression of 

In(M) = -21.49%C ,2*2*f5/3*)-1/3 (2 3) 


where C,* is the only parameté2r. 
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III. EXPERIMENTAL PROCEDURE 


——e Pee ae ee ee ee a ee oe a == 


The experiment is perform2i using two different laser 
sources, helium-n20n (He-Ns) and gallium-arsenid?e (Ga-As). 
Figure 3.1 shows a block diagram of the experimental set-up 
used for both lasers. 

The measurement 9£f C,? along the optical path is made by 


-~ + 5 ‘ Ca) 


feeds <¢ VidicoOn and t3.1escope 25 tA 
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Mimene Dasement of Spanagel Hall. Atmospheric turbulence is 
preduced by nine overhead hot air ventilators. The optical 
Semepment Similar ¢o that described in reference 3, consists 


of a 6 inch diameter Cassegrain telesc 
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focal length [Ref. 3]. 

Meeng netral dansity absorption filters to attenuate 
the intensity of the laser bean, the telescop= is 
Meeimanated at its liput apertur2. The light image from the 
telescope is split by a beam splitter with one beam sanz to 
fmmemerealcOn and the other to a slit scanner for comparison 
ert the two systems. Ties. 25 Then Eaanstormed [oO an 
anaiog signal and recorded on a Panasonic(NV-9240) video 
tape recorder for processing. THe Vata son thst slic 


scanning system ar2 recorded on 2 precision instrument “ape 


if 
1) 
| 


recorder for analysis. At EIcording the video signal, 
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memeysoS iS accomplished using a T¢sktronix 468 Digital 


J 
Seepage Oscilloscope and a Healstt Packard 9825 Computing 
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System. The linearity nasonic tape recordsr was 
demonstrated in Crager's thesis. This is also confirmed by 
observing the images as real time and recorded 
@eepla ys [| Ref. 3]. 


The 468 oscilloscope is set-up in the following manner. 


Peeanpuc Signeil to channel A is located on the scope using 


jene non-storage mois. Pie sagmal 21s thon .displeved wsing 
mmemA LNTENSITY switth of tn horizontal display. Peas 


Mimeens risa zone is used to position th B sweep (delayed) 


(i. 


so the desired location within the A sweep interval to 
obtain an expanded view of 2a waveform for 2xamination. Once 
Sem wWavVelOrm is seatered on the scope, RPE’ “Mer LZoule | 
Meferay -=S5 Switched *> B DLY'D. Lio soea One to f2Cn lt ta ve 
Pema: gital storage Slrcuitry time base by using the satting 
Mieeenc B TIME/DIV switch. 

The waverorm is 2xpanded in time by decreasing the A 
TeieyDEY switch setting and moving the waveform back to the 
Middle of the scope using the delay tims positicn control 
dial. When a representative vavetorm is obtained, eedicn 5 


femey7WiY switch is then used in conjunction with the delay 
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mfeesingl= central naxkimum of the waveform is centez|ed and 
one horizontal sweep is display3id. 

Once this waveforn is satisfactorily obtained, the AVG 
storage modé is sealiscted. Th 468 will average the input 
Signal for a selectad number of sweeps and display «he 
accumulated waveforn. All data for this thesis us2 32 
Sweeps for éach average. Hien Ooms 1 DOWN Ecady “= Cltranster 
Gata when interrogated by the HP 9825. 

ineemcOMpPUtTEL 15 he coatrolier for ali intertacing 
operations Wilt pie OMY 42 miiinal amount of operator 
moeeerac clon. Diem eho Oeea tee trace Ss Maauly “cOoten suss 
that the equipment is properly set up and to select if data 
Beento b2 plotted. DEGLtUZac2On Of ce analog ¢ Signa cs 
meeempizshed by the [2 ktronix isos Processed waveforn data 
are transferred fron the misroprocessor mémory <to the 
Beeregs Display RAM [Ref. 8]. 

Mae controlling scogram of the HP 9825 interrogates the 
458 via the IEEE 488 interfacs> bus. The 468 receives the 
Gata request from the controller and sends the waveform 
message, both preamble and dat:. The waveform message is 


Seered in the calcuiator memory ficr further processing. 


‘ 





When the message is complatea, the 468 concludes with an ena 
Sens couction termiretor and the controller takes control 
of the bus again. 

Bata processing 2f the digitized signal begins with the 
Gata being stored on magnetic tape for further use by the 
computer. The data are stored, processed and plotted by th 


Main program. Subroutines ar2 stalled as necessary for their 


specific uses. Th2 two signal waveforms are czecorded and 
processed. First, the signal from a laser beam inciden+ at 
the aperture of the t2lesctope is recorded. Figures 3.2 and 


3.3 show these recoried iata waich are rafarred to as thks 
Sseurce data. Seconi, EOaMCT LOT ati On, the signal from a 
Me@mmeoae md With 2 adiftractidnh grating in solace 32: ths 


aperture of the telescope is re 


(2 
© 
r4 
Qs 
ch) 
Qu 
r 
tJ 
ra 
(D 
(D 
= 
ct 
wy 
rf 
Mm 
" 
mS 
in 
ct 
(p 
| 
be 
OF) 


Beeeeobrated by uSing a grating that consists of closely 
Spaced vertical bars in front of the telescope. This 
produces a diffractiosa pattern in the image plane and allows 
memmeen= Caliculation »9f the scale factor. These recorded 
data are referred to as the scale data and are plotted in 
Peouees 3.4 and 3.5. 

BiGeng the recorilng process, the signal from the laser 


beam is recorded on video tape while the signal is obs 
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often very sensitive to minor movements 9f either. U pon 
@emetetion of recording the data on video tape, the 
experiment is ccncerned with th2 compilation and analysis of 


gata. Mme LLEsSt ieee to be determined is the scale farztor. 


iD 
roe 7) 
~@ 
an 
cl 
=; 
fee 
f- 
t— 
1} 
D 
= | 
su 
i? 
DP. 


Once the scale factor has been i2termin 
Gem@erant throughout the calculations, unless tne Telascope 
Meet length is changed. 

Mrenscaie factor is calcul2ted in tne following manner. 
The mumber oF points between the peaks in the plot of the 


Mee@ien dasa (Gijtfraction grating in solace in front of hs 


telescope averture) is measureai. wnen the distance between 


peaks, the spacing between the bars in the grating, and the 
wavelength cf the iaser are known, the scale factor can be 


calculated from the relation 


Sint ehet a) = Ara (3. 1) 


Sin (theta) ~ theta (352) 


(Small angl2 approximation) 


sft 


(cheta)/ndp C313) 


st A/ (d*ndp) (3.4) 
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scale factor in radians per point 


7) 
th 
i 


A = wavelength in n2aters 


d = spacing between lines on jrating in meters 
ndp = number of points between santral maxinun 
peeomericst Order diriractcisn peak 


Thas standardizes the data for the abscissa of the plots in 
mts enS p2= point. 
The program now takes the transferred data and computes 


Muppereene Spread functidn as shown in FPagures 3.5 and 3.7. 


Seeet en ~egqration using squation (el) Wes: “OD 2Ge09 G. «ane 
Spread function. Nee eho robo er ~transrorm Of che Line 
Spread function is calculated. These curves are plotted in 


hegures 3.10 and 3.11. 

jimenmd=tiraction of the optiss is comput2ad and plotted as 
shown in Figures 3.12 and 3.13. Following the same mathod 
Memeeeaore, the line spread function of the optics function 
Bemealculated and plotted in Figures 3.14 and 3.15. The 
program now takes the Fourl2er transforn of the optics 


mumercaon LSF. Thes=2 data 2re plotted in Figures 3.16 and 
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Ball. Pro Pelllract SoenctOrnNs Swe che system anid cotses are 
Mevecded point by point t> yieli the MTF of the atmosphere, 
meeesonowr in Figur2s 3.18 ani 3.19. Pie be» Of ene 


atmesphere will be usad t> calsuiate the long-term value of 


Smeets calculated by linear regression of In(MTF) versus 
#5 /3, Ere SLOpe Of Ehe regrssston is proportional to C,°. 
Two Simuitaneous 2quations ara solved using Cramer's rule 


and the matrix ROM of the HP 9825. The equations used are 
a*X + bD*¥T = Y (355) 
a*X¥2 +4 bH*eX = Y*x (3.6) 

where 

a= slope of the curve 

pe) intercept 

X = Summation of (I*sf) **(5/73) 


X**2 = Summation of (I *st) ** (10/73) 


I = point number (total number = 256) 
Y = Summation of natucal logarithm of MTF 
Y*¥*X = summation of products of two values 


Zo 





Mees Cbtained from the apove inrormation using equation 
Wes). This will yield 
Cyt = af (-21. 49*Z *£** (5/3) * ** (- 173) ) (S74 


jme Program now goes to the prediction phase after 
Gaiculating C,# f5c the atmosphere. Tf a Gaussian 
@eerteibution for an input sours3s and a value of the standard 
Meweation are assunad or “known, ene, Sesulitang Wsecurce 
Mire... On can be calculated. A DaLOe “Cl etn s scomputed source 


MeesrOwn in Figure 3.20. It caa be seen from this olot that 


Seg 


the half-width at half-maxinun is approximately 
MeelOradians. AS DSF OuS,) ane Line «Ssopead TUunCceiOn OL che 
Meee ee SOUTSCE iS calculated a234 plotted. This is shown in 
Megi2e 3.21. Phe  POUrteSr. “LreasroOrm. “32 Jthos | “dacs ies 
Merculated and plott2i in Figure 3.22. 


mre program now multipliss the Fourier transform of th 


() 


computed source with the transform of the system (includin 

the atmosphere) ani plots the result in Figures 3.23 and 
B24. Peots Of the Lnverse FSurier transform are shown in 
Mmegueres 3.25 and 3.25. Next, the Abel transform is comouted 
Meme lotted in Figur2s 3.27 ani 3.28. The Abel transiorn, 


as described in Chapter II, transforms a one-dimensional 
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Mere Spread rinction in 
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Boeus R is calculatei and plotted as shown in Figures 3.29 
ema 3.30. Diss LoaeNe GEastlon Of power that one would 


Gd value 


(ip 


expect to be incident on a tarjy2t using the measur 
of etmospheric turbulince as an input parameter. A detailed 
Pimaevys.S Of the computer progran is the subdject of the next 


emapter. 
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SOR Lees 


LANGUaGeE 


— oor cee ae ae ee «= ow 


fie Drogsam pPresei ted hetS 25 Written in the aigorithmnic 
danguage 2s outlinel by Grahan [Ref. 9}. inc vic Ge] 1ded 
that the language in use here be machine independent. 
Comments appear in parentheses. 
Algorithm Laser 
ferpu~ which laser is used.) 
Input Laser_davelength 
ESTRING <== Type_dsf Laser 
er put 1 MileMbaretellar tasst Source data --ere recorded: 
mre fOllowing lines letermine if the laser source dati are 
already reccrded. Pees, SOULS HS bda-e ers Got TecOrd ed “tier 
Mme DrOogram goes tds the transfer subroutine to 3llow 
MemeerctecangG With che 468 to seflect a data transfer. The 


parameters passed to the Plot subroutine h2 


Seogniiicance: the paranets 


Tie 
— 


determines 


ve the following 


2ho HOt. Zones | 


Sraex, Cor example 512 points or the number of microradians. 
The second parameter tells the computer which file the data 
are recorded on so ta2 data can be loaded into memory. The 
final parameter determines the labeling of the horizontal 


- 





—m «2 ws 


mecroradizn.) 
If Laser _Source_not_recorisi then 
Call Transfer (Source_Fils) 
Cater TLOttol2 ,sousce, File, 0) 
Bea if 
Meeewe 2f£ the particular laser scale data é@ére recorded. 
fewer deta are recorded with che diffraction grating in 
place. Tf the scale data n12ed t0 be recorded, the 
subroutine Transfer is called and the scale data are stored 
On a separate file.) 
ieee bases Seale Noweresorded chen 
ea 2 i tenstan (Seale & 
Car) Perel? oes ls fils, 0) 
md. Lf 
mut Scale Fastor 
fee2culation of the scale factor is detailed in Chapter [II. 
Comparison or +the time base of the recorded data is 
necessary to ensure that the data recorded have <zhe same 
time scale on the oscilloscop:. If the time scales are 
different the program stops ani displays an error message. 


A file is loaded, *h2 waveforn oreamble is searched for the 


wo 





recorded files: source, S24 )57% and 
Colliimated data refer +> the same laser 


maen=rolled conditions. The laser bean 


Pewensure consistency of invuts.) 


Input Laser_source 
B<== pos(Laser_ Sdurce, *R:') 


R1<s= val('Be2',B+4*) 


PoenaiNG<== Lassr source ('Bt3','s+u") 
Input(Laser_Scal3} 

B<== pos(Laser_Scale,'R:°*) 

Rew-— Val('St2! ,3+4') 

B<== pos(Laser_Stale,'IT:') 
BSTRING<== Laser_Scala('B+t3','B+ed') 
Input(Collimated_ source) 

Ba-= pos(Collimatad Source,'R:") 
R3<== val('B+2!','8+4) 


B<== pos(Collinat2d_Source,'IT:") 


2)! 


say) nen—litveeritaes of th] video recorder, 


BRUnITS" CES UM XINGa 5) Vee th shop eleyet= Settings orale 


Peec2 cular data ar2 stored. hes 2S "done. TOr.e@ack Cz 


sollimated 


is projected 


under 


ina 


@omlcamator and recoried on tape. This data could be used to 
electronic 


Semepment, Or recelLvlag optics. Data sl2anents are compared 





Pooks ntNG = Sola LNG = CSTRING 


Input KkDOATA(Fil2_Number) 


Rae Le 

Base 
Oie-puT "ERRIR IN Prete Scales 
Saw eps, 

ema Lf 


YDATA (1) <== RDATA (1) 


Ema Dod 


Mir=s section takes the data 


the 468 and 


meensftsrred =rom 


Souuecs ( R43) ' 


Bey) 


heh 


from the waveform message as 


removes any DC background from 


an image as well as subtractiny the minimum value in order 


to "zero" the data. This makes 


mee pOth calculation 
and greater than 463 


This represents the 


ang plOurL Ad |pus poses. 
are sumned and stored 


(Wetgsce Of... e. CULV ec. 


V=Ouas chs Sah sing Solint 


Data up <9 517 
20° 3D. 


Phies 


process removes the DC background from an imagé.) 


Do for I = 1 to 512 





aA<== aAtYDATA(T) 
feo ‘<2 Sd ea Sh 

B<== YDATA(I) + B 
ime .f 
If I > 463 then 

B<== YDATA(I) + B 
eat £ 


ma Do 


| 


Meutput Cor debuggi2gj) 
Ouzout A,B 
(Average of the sum 9f data vaiues in the "wings".) 


D<== B/100 


(D 


Menas is the total Sim of all data less the “wings” value.) 
A¢== A-512*D 
Meet dat2 are summeil and storai in C. Each data 
Bepcompared with th previous total dividsi by 2. 
Mres Will find the index value, I, for which 
the maximum value exists.) 
Do for I= 1t)9 312 
YDATA (I) <== YDATA(I)-D 
c<Cs= YDATA(I) + 2 


Piece 2/72 then 


ty 
\O 





th 


Bee 1 
moc Do 
(Output for debuggiag.) 
Sueput C,J 
(The curve is shift2di so that th maximum value starts 
femmes Origin of ths horizort2l axis. This is called 
the point spread funzttion and is recordei as RDATA.) 
Bomeor K = 1 €)9 236 
Tenn ACh) <= — 
YDATA( (J-K+t 1) moad512+1) +YDATA ((KtJ-1) moas512+1) 
End D> 
Mentor T= 257 to 512 
seit os ) 
mma Do 
BieetOor 1 = 1 t5 256 
RD ATA (I) <== IDATA (TI) 
moa Do 
(Determines if tabular or plotted data ars desired.) 
If PSF_OUTPUT_wanted then 
Call Table (RDATA) 
Ema Lf 
Prt PSPePLOT wanted then 


Sa ie LOt (32, RDALA,. 1) 


a0 





tty 


Bod & 
Meetcttiaces the line spread fuastion from the poinc= 
spread function.) 

Gatl LSF 
(Determines if tabuiar or plotted data are desired.) 

mE Lor OUTPUT wanted then 

Call Table (RD ATA) 
Besse Lf LSF PLOT _wanted thea 
Geom” Plores2; ROATA 1) 

Emad Tf 
Meerculaticn of Fourier [ransformn of LSF.) 

Hepuct RDATA 
(The maximum value of the array IDATA is at I=1. Th 
Mmempu- back into array RDATA 2nd folded over so that 
meer A (257) <5 IDATA(512) is th2 mirror inage of 
meee (1) to IDATA(255). Then this array is recorded 
as RDATA and the Fourier transform calculated.) 

Peetor I = 257 td) 512 

ALA 1)<== 
Emad Do 
Do for I = 1 to 256 
RDATA (514-I) <== RDATA(T) 


end DO 





(Determines if tabular or plotted data ars desired.) 
If FXFORM_OF_LSF_JUTPIT wanted then 
Call Table (RD ATA) 
Else If FXFORM_OF _LSF_PLOT wanted then 


Call Plot (32, RDALA,2) 


rh, 


Emo. 
(Dia of Cbscur/Di2 of Joj Lans in meters, RATIO=.064/. 164) 
(This is the measurement of the ratic of 
diameters of Sugeneanecen Eo) co ese opeS.) 

tal= Scale Of data in microradians=1.60) 

(Diameter of teleszsope Obj Lens in metars,OBJECT= 
Q.164 neters) 

ime Laser Source = He-Ne then 


Wavelength<== 6.328%*2-/7 


ty 
}- 
©) 
(D 


Wavelength<== 9.95*e-7 
moe Lf 
(This begins the salsulation o&€ the Fourisarc transform 
of the optics uSing the Pay Seimeccren..at! of =he 


data regarding the optics are recorded on a separate file.) 


thy 


(These next four statements fic the constants in «the 


a2 





meagument. Of the Optics function.) 
beli<=—=) Scale Fackor 
D<== RATIO**2 
H<== 1-D 
Z<== PI*Scale Fastor*Ibject*1*e-6/Wavelengin 
Ber fOr L£ = 1 td. 256 
Y<= Z*(I-1) 
me Y > 30 Enen 


RDATA (I) <== 


iJ 


+ 
7) 
(D 


RDATA (I) <== (AIRY(Y) -D¥AIRY (Y*¥RATIO) ) #¥2/4 **2 


Siva © f 
K<== [I 
Eada Do 


(Record RDATA on s2parate fils) 
Mtmeanmines if a2 plot of the optics function is desir2i.) 
Iit OPTICS FUNCTIIN_ PLOT wanted then 
Gail Piet (32, RDAT A, 1} 
ead If 
feelculates the line spread function of che optics 
mamSct ion. } 
Sarl LSF 


MPe=ermines if a plist of the LSF of the optics function is 


6 3 





desired.) 
Merlot OF Est GPTICS wantei then 
Call Plot(32, RDATA, 1) 
pad If 


Do for I = 257 ¢> 512 


RDATA (I) <== 
Emc Do 
Meee re £OLGEC SVSt =2C63in 2m tae Feurlis: *transicen 


caiculated.) 
Hemet Oor J = 1 to 236 


RDATA (514-1) <== RDATA (L) 


id 


ea Dd 
RDATA (257) <== RDATA (256) 
@eeGOrd xADAT A) 
Beeuriet troansrorm Of the LS& of Ootics is calculata3i.) 
Sent ¢X FORM 
(eetermines if plist of Fourier transform of Optics 
function is desir2i.) 
me PLOT OF_FXFORM OF _JPTICS wanted then 
Gal  PLOt (2554 8 DATA, 2} 
ead Lf 


ee culation of th 


(p 


Gigt@a=f. of CWO FOUTLEE transforms.) 


(This section takes the Fourier transform of the systan 





Mamet vVicsS =“u, polnt GY polac, by the Fouriss 
of the optics. This yields ta12 modulation trans 
Function (HTF) of the atmospaere.) 
Input RDATA_OF_ FXFORM_OF_SOURCE 
BO tor T= 25/ t> 512 
RATA Ci) <== 0 
End Do 
BoecOr © = 1 to 236 
RDATA (514-1) <== RDATA(TI) 
End Do 
Input RDATA_OF_FXFORM_OF_JPTICS 
BewtOr IT = 1 to 236 
RD ATA (I) <== IDATA(LI)/RDATA (TZ) 


Bhd Do 


Meererminss if a pidt of the fF is desired.) 


iweee bot OF Mar OF ATMOS wanted then 


Serie OU Ic, ADA LA 1) 


MepeculatiOn Of Cyc by curve fitting. 

mies Procedure is iescribed in Chapter III.) 
(All "R" variables used in this calculation are 
Meicecalized to 0.) 


R1V7<==kR 186 ==R19C== R20C HS=R21C = =R22<KS=R23 C==R24NK 


Po 7 
Eee 


22.97 5c os 





Meee ysis cr data shows that computer «ine 
Beyond 96 points.) 
MOmeEOr I = ts 95 

fee ~ when summed is ths (f(1))19/3. This is the 

total of the angulac spatial frequencies squared.) 
Ris<=— ((1-1) Scale Fastor) ** (1073) 
R17<== R17 + R18 

feo - when summe? is “ne (Z(2;) 5/3. This is 

mime total of the angular spatial frequencies.) 
R1%K== ((I-1) *Scale Fasttor) ** (5/3) 


R20<== R19 + R20 


Meee 25 the =ctai number of ooints.) 


Do Whiie RDATA(I) > 9 
(R22 - when summed 15 the tot32l of the product of the 
Mteural Logarithm of each point with R19.) 

R22<== Ln(RDATA(TI) ) *R219 

R23<== R22 + R23 
(k24 - when summed is the toral of the natural logarithn 
Sr C€ach point.) 

R24<== | n(RDATA(LI) ) 

R25<== R24 + R25 


End Do 


Oy 
O*V 





(Outout matrix values for debuyging.) 

Output KRIW7, R20, Rk21,R823, R29 

(Assign value to matrix. 
This section takes the above calculated values and 
sets up a matrix. Phe natrix 2qguation soived is 
Ee=— A-?*C,) 

Peres tes Veet lS Ss A (2/2) 401241) go (251) -) 

A (1,1) <== B20 

A(1,2) <== R21 

A (2,1) <== R17 

A (2,2) <== R20 

INVERT MATRIX_A 

C(1,1) <== R25 

C (2,1) <==R23 

B<== A*X¥C 
(Output matrix values for debugging.) 

Smweput B{(1,1),B(1, 2) 


fetreulation Of C,° using squatica (2. 3).) 


R22<== B(1,1)/(721.49*Range*®R20* (Wavelength** (-. 3333 3))} 


SucDut 'CNSQ=', R22 
(The program now starts the prediction phase. It 


calculates a sSourt2> assumed => have a Saussiaén 


5) Hi 





mee program uses the calculated vaiue of C,* and predicts 
the power incident on the target.) 
(Computed Source Pattern: A=A)*exp (-x**2/2*Signa**2) 
(This is an "arbitrary" Gaussian source pattern 
with a standard deviation for Sigma set squall tc 2.) 
Me<== A ZERO 
C<== Sigma*Scale Factor 
Peetor I = 1t> 2536 
Poo Bibi) soc ale Pract o> 


G<== PR*2/ (2*C¥*2) 


RDATA (1) <== 0 

Bod DO 

RD ATA (1) <==A) *eExD (-G) 
End Do 
(Record RDATA. This records the computed source.) 
(Determines if plot of computed source is desired.) 
meee LOT Of COMPUITED_SOURCE wanted then 

Sa di ePlOc(S2 2k DALA, 0) 
ema Lf 


(Calculates the line spread function of the computed 





Securce.) 


aut LS 


Me@ececrmines if a plot of the LSF of the computed scurce 


is desired.) 


If PLOT_OF_LSF_OF_COMP_SOURZE 


Call Plot (32, RDATA,0) 


DO LOr 


RDATA (I)<== 9 


L =)25 7a. 512 


RD ATA (S14-1L) <== RDATA(LI) 


ena Do 


RDATA(257<== RDATA (256) 


eat FX FORM 


1 then 


M@eerccmines if a plot of the Fourier transform of the 


LSF of the computed source is jesired.) 


If PLOT_OF_FXFORM_OF_ COMPUTED SOURCE wanted then 


Call Plot(255 ,RDATA,2} 


ena If 


(Caiculation cf th:a 


product sf two Four 


3s 





Weear = pe MIF CE she atmosphers. 
Bo fOr [ = 1 t0 236 
RDATA (1) <== 
RDATA (I) COMPUTED SOURTE*RDATA (I) _FXFM_OF_MTF 
End Do 
Perodiuct of FXFORMS of Source*Atmosphere*0 ptics.) 
ADATA (1)<== 
Mien ') PRODUCT JF 2 FXEFUS*RDATA (1) FXFM OF _OPTICS 
Bortor [= 2 t9 256 
RD ATA (I) <== 
Se Ame ne Pee Or ye 1b XP MSeRUATA(.) FAFU_OF. OPTICS 
ee ee. en 
end Do 
RDATA (257) <== RDAIA(256) 
WeeeweDove data predicts what the Fourier transform 
Seeene entire system is, using the calculated 
value of C,* for the atmospheric «urbulence.) 
(Determines if a plot of the result of Fourier transforcn 
is desired.) 
If PLOT_OF_FXFORM_ PRODUCTS wanted then 
Saet - hOt (32, RDATA , 2) 


end if 
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Gye =Se FOUEZCS THaNSEOEMN Jives taraqet LSF.) 


Gave TUVER SE FAPIR | 


(Determines if a plot of the iaverse Fourier transform 


is desired.) 
If PLOT _OF_INVERSE_FXFROM wanted then 
call Plot(32, RDATA, 9) 


mmo if 


two dimensional PSF by Abel transform.) 
Call ABEL 
(Determines if a plot of the Abd2l transfora is 
If PLOT _OF_ABEL_XFORM_wanted then 
Ga it” Paset (32, RD ATA, 0} 


End If 


feoentla tes the fraction of odwer insides circle of radius 


Fem inhas predicts the fraction of power that 
mimeadent on target. } 
RDATA(1)<== 0.25*PI*RDATA (1) 
Do for I = to 236 
RDATA (I) <== 2*PI*RDATA(I) +RDATA (I-1) 
End Do 
(Determines if a plot of the fraction of power 


If PLOT _OF_ POWER wanted thea 
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Cae Ploc( 32, RPAL A,)0) 
Ema LL 


Enc Laser 


(The following subroutines aire used sl ganas 


Laser.) 


D 
= 
G 
=) 
Oo 
(w 
ie 


SaoLoutine Transfar(Fil 


he Algorithno 


(The transfer subroutine gets caw data, both preamble 


V2 
jv 
iQ 
(Dp 
UI 
Mm 
rie | 
i. 
ou 
me 
re 
Lo By 
(Dp 


and data from the waveform mes 
Mee. ©£t processes the wessag= by finding the a 


value of the array then subtracting-this value 


2nimum 


ra on 


Seen SLiemert in the array. This "zeros" the array.) 


Dieout "SHSULS EQVUipmSnt Set up properly’ 


Output "Continue #hen reaiy' 


Do While E = -1 


("DATA" is the I/0 buffer where the data fron 


Sent. The status of the buffer is read while da«a are 


transferred. When th2 transf2r 1s complet3a, <= 
is cleared.) 
Data<== Trans ferred_Data 
ieee cdelhanster Compisce chen 
E<== 0 
(2) 16h ely 76 3d 


[Le 


the 468 





Clear _Interrace 
N<== 1 
(This is where healer format stops and data beqin.) 
A<== pos(Data,' 4%") 
Do for I = At2 to 687 by 16 
Do for J = 1 to 16 


LT&£ J+I>588 then 


Oo 
Ly 
«nN 
it 
it 
QW) 


Ne== Ne 1- 
EenGat © 
RJ<== num(Data(I+J) ) 
(The numerical values of each szlement are stozted in Y¥DA 
YDATA (N) <== RJ 


N<== N+ 


E<== mininum(YDATA) 
Demcor L = 1 to 512 
YDATA (I) <== YDATA(I) -z 
RDATA (I) <== YDAT A(T) 
End Do 


Store_RDATA (File_Number) 


3 





Re=v0rr 


Pea Transt er 


(This subroutine calsulates th: 


for the previously r2acoried poi 


by using equation 


Subroutine_LSF 


(Depending on which file number 


feon the calling subroutine, 


data are loaded.) 


Ir File Number 


(1. 1) <) 


Input Scale Jata 


BD) (MESSE Ce NAVE 


Else 
Input RDATA 
acd «If 
Beetor £. = 3} to 512 
io Maen (030) Qe 
mma DD 
Do for I= 1¢t95 256 
para tt) <== 
ena Do 


(Plots have shown that comput2ar time is 


valuable information is gained 


2 DATA (I) 


2ither scale 
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bame Speead Funcrv=ion 


nt spread function 


has been pass¢d down 


Je oouL Ce 


scale Data then 


Lost and no 


Devyona avo. 24° DOLNts.) 





Doster f a E24 
J<== 1 
Output Lf 
Q<== IDATA (I) 
Do While Rk < 24 


R<== SQRT(L*I+5* J) 


Geet t= eres Ul ona ss (20s (2) ) Fi reccicn (sp =2 (242 (2) ea 


J<== Jr] 
End Do 
IDATA (I) <== Q 


mad Do 

Be Lor FT = 1t5 24 
RD ATA (I) <== IDATA (T) 

End Do 

Mmeeraie Number = Scale_Data then 
Store _ Scale Data 

Else 
Store RDATA 

Baad If 

Return 


mad LSF 





feris Subroutine calculates th2 Fourier transforn 


of the given data using the Cosley-Tukey Algoricthn. 


Tf the inverse statement is true, then ths 
Fourier transform is calculatei.) 
Subroutine FXFORM 

Set_Radian_Mode 


(2**9=512 which is the numbec of points.) 


ber ror If = 1 #3 5312 


ID ATA (I) <== 3} 


(24 


nd Do 
If File Number = Scale Data then 


in Plceescal ce.) ate 


Input RDATA 
eg LE 
T<== PI /2** (N-1) 
DOmrO r= 0 tO et {4— 1) = % 
fee2S a2 bit LnversSidn subroutine.) 
Cav ei ( J, ?,N-1) 
@<== 655 (P*T) 
Iz INVERSE FXFORM = True chen 
Fl3j7<== 1 


76 


“nverse 


: 





rij/7<== 0 
Bud 25 
P<== sin (P*I) *(1-2*F1gq7) 
DO DIEh = Zee BIt) cto ZeROeI+ RO 
RPG =—— RR DATA (2) 
R2¢== RDATA(I+R0O) 
AG] TOATA Ge) 
R4UC== IDATA(LI+#RO) 
RDATA (I) <== 21+R2*C+R4Y*D 
IDATA(I) <== 23+R4*C-R2*P 
RDA TA (I+RO) <== RI-R2*C-RY*P 
DOP A(ET RO) C== ROR R44 CHR 2 SP 
znd DO 
ena Do 
OED tr Mh 
End Do 


ums Section is for re-ori 


nu 


Seed 


—_—_ »_ as 


NG ene DL OOK <) 


0 to 2* *N-1 


Sa aL (i, dy ) 


DO 


While I-J < or = ) 
Do Whils I # J 
PCZ== RDATA(IT+1) /SQRT(2**N) 


ia 





<== DADA Vy SSORT (OMEN) 
RDATA (I +1) <== RDATA (J+1) 
IDA TA (J+1) <== IDATA (J+1) 
RDATA (J +1) <== P 
IDATA (J+1) <== Z 
Ee DS 
RDATA (I+ 1) <== RDATA (I+1) /SQRI (2**N) 
IDATA ([+ 1) <== IDATA (I+1) /SQRTI (2**N) 
Ena Do 
If File Number = Scales _Data then 
Store Scale Data(File Number) 
Hise 


S=0r 2 Sara cA Ce 


i- 
(bp 
L, 
bs 
hase | 
= 
(D 
ty 
Sew 


fa sofel, 92846 

moa D> 

Set_Degrees_Mode 

Return 
sna FXFORM 
(This subroutine takes che binary number P1 containing 
Pemoets and inverts it end for end, ¢.g. 3910111 becomas 
m1010.) 
Seo eoutine BI (P1,P2, ? 3) 


P2<== 
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POowteam. = iI cS P3 
PYU<== P4/2 
Peae=] 24D? 


Meee rect ont ou) & 0 Enan 


P2<== P2 + 1 
Bed [2 
PY<== Int(P4) 
Boa D> 
Return 
Bnd BI 


Weees LUNnCtHen calculates AzEy(x) =2*J1(x) /x whete J1 {x} 
es the Bessel function of order cne.) 
munca: On ALRY(P? 1) 
fer 1 < 0 then 
Cute. “ERROR -AanGueew tl Less PHAN J$ 


Sc Cp 


Return RY 
End If 
RS<== 


yo 





Me Pl > 15. then 


R6<== 90 + PI /2 


R12== 1.4*P1 + 60/P1 


Emad rf 


meeP) < 5 then 


R6<s= 20 + 10 ¥*P1-P1**2/3 


R12¢== 6 + P71 


Else 


R6<== 20 + 19 *P1-P1**2/3 


R12¢== 1.4*P1 + 50/321 


moe LNW (2 129 ne (3+ 174) ) 


Meee teo es Oo Oy 3 


R8<== 1¥*0e-28 


R13<==RU<== 0 


If M/2 = int(M/2) than 


FligiQ<== 0 
Else 

Plol0<== 1 
End ££ 


Do for J = 1 to 4-2 


R1S<s= 2*(M-J) *RB8/P1-R13 


R13<== 28 
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Ti M-J-2 = 0 then 
R4¥C€s= R15 
End If 


If Flgi0 = 9) then 


Ble Oc=-= 4 
Else 

Fljy10<== 9 
Ged 


RIGC== 2 144+2*RB*F1LI10 
End Do 
R1X== 2*28/21-Rk13 
RI4<== KRWUtR15 
R4<==R4U/R14 
Tf (abs (R4¥-235)-abs(R&¥*1*2-6)) < of = O «hen 
R4<== 2*R4/P 1 
Eyl. L & 
Retucn RY 
R5<== RY 
End Do 
weeput ‘ACCURACY NOT JBTAINED' 
Return R4& 


End AIRY 
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frre subroutine Abel takes a 312-dimensional 
line spread function and calculates a two 
dimensional point spread function.) 
Subroutine ABEL (Fil2_ Number) 
N<== RDATA (1) 
RDATA(1)<== 1.4*2DATA(1)-1.3*RDATA (2) +. 4*RDATA (2) 
DemerOr -< = 2 t5 63 


M<== RDATA (I) 


RDATA(1I)<== . G*N+. 24-. 5=RDATA (7+1) 
N<== ff 

End Do 

WOetOor f= 145 5% 


ee bam S—==— 95 DADA (2) 7 (2-S ORT ( (1+. 1) = *2—LF Tt) ) 
DOsroG ws = 1 <o.54 
R DATA (I) <== 
Poe eA) +R DATA (J) /SORT ( (J+. 1) **2-1*T) 

End Do 
RDATA (I) <== RDATA (I) SPL 
Output I 

End Do 

Return 


End ABEL 
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Mons = SUbECUTIne Prints out 1at3 in tabular forn.) 
Subroutine TABLE(Fil? _Nunber) 
Input RDATA (File_N umber) 
Boeros £ = 1 t> 32 
Dor tas. J = 1. tocls 


Jutput RDATA( 16 (I-1) #J 


Output RDATA (16*I) 


End TABLE 


Memes Subroutire plsts the desired data.) 
(Each time poets call ea, tls Subroutine 
Peeks the particular qraph and labels ic. 
Mm@emaOlTizOntal axis is labelid according t9 the first 
and last parameters passed by the calling subroutine.) 
Subroutine PLOT(P1,Fil? Number, P3) 

Input RDATA(File_Number 

A<== 0 

B<== PI 

C<== minimum(RDATA) 

D<== maximum(RDATA) 


XMIN<== -. 1* (B-A) 


3/6, 





KMAX<== B+. 05* (3-4) 
ViN<== C-. 1* (D-<c) 


YMAX<== D+.05*(D-2) 


E<== 8B 
F<== 10 
(Using the value of Pl, eyo en OE?ZON TS. TECrement «fo 


plotting and labeling is jietercmined.) 


3<== 8 
Else 

ee= 
mma Lf 
Eeot B,C, 1 


Merton i nes and plases tlic marks on the horizontal axis.) 
Deoetor f= 8 t> 0 by -G 
Elec <b Al 7 CtA, C, Z 
PebOct (b= Ayy ot A,ct(D=- 5) 7 150, 2 
Piotr Gp—- 4) 7 p+, C, Z 


End Do 
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Plot A,I*(D-2)/F+4C,2 
Plot A+ (B-A) /150,I *(D-C) /Ft=, 2 


Plot A,I*(D-C) /Ftz,2 


tr 


na Do 


@aaracter Size 1.2,1,./7,0 


(@enas numerically labels the vertical axis.) 


Beeror f= fF to J by =-1 
Baler eh Ov ott Ob - Aja. (D- Cy) JP tC 


Label I/F 


If P3 = 0 then 


ASTRING<== "DATA POINTS" 


re 


ASTRING<== 'LINES/MICROIRADIAN' 


L< == 1/(2*1) 


Ena If 


(This numerically labels ths horizonta 


De zor Lt = ,0et> E by 3S 
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PERODmAtT (hen. U2 oO) (o-Ayn, o> .6 025% (D-c) ,1 
Label I*L 
End Do 
Tf P3 # O then 
L<== Scale _ Factor 
Ema If 
f1mrs labels the horizontal axis.) 
Plot .4*(B-A)+A,.)5* (D-C) +#C, 1 
Label ASTRING 
Plot -.07* (B-A) +a, .3*(D-C) +c, 1 
mriaeact er Size = 1.2,1,.7,90 
ASTRING <== "NORMALIZED INTENSITY' 


Label ASTRING 


Plot I, maximum(RDATA(I)),1 
Bo for IT = 1&9 PY 
Plot I-1,RDATA(I 
Ena D> 
If P3 = 0 then 
Input PLOT _LA BEL 


ISTRING<== PLOT_LABEL 
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“J 


£ 


Ze) oa 


t 


Plot «Ob, .9*D, 1 

Eaocl * £LOT OF! 

Pao. 2.695, .57*D, 1 

Label ISTRING 

Poe 2648, .S4*D, 1 

Label ESTRING, ' Laser! 


Return 


Fra PLOT 
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Ve weeowenos TON 


ivessWwOLK LepOoOrssi in this thesis supports the aodel 
predict2d and measur2i in Crittenden, ana ewethers.» for ths 
long exposure case [R2f. 2]. ite vids cor, 1m ceplacing the 
mechanical slit scanning syst2m, shows no degrading of thse 


Seaman data. It also demonstrates that a good approximation 


a 
sv 


FOr the point spreai function aav be made Dy recording 
Seqgte [TV line through a lasér spot. This Line is then used 


to calculate the one-iimensional line spread function. The 


jinearity of the vidso tape recorder is s22n in the rasults 


Semone MIF and C,@. Mie soc DoculicScop> 2S the workhorse 
MOnecac entire syste. mono he-Care..y 4dr solays, “Stores, 


The overall syst22 does not have the capability of the 
system described by Cragerc. However, the relative 
Sailipuzec.ty cr the structured programming t¢echnique coupled 
with a lower equipnent cost ismonstrates tnat comparable 
measurements and data evaluation using this equipment can be 
made. Further investigation shoayld include use of «he Data 
Precision 6000 digit321 waveform analyzar thereby allowing 


data <«o9 be sampled and analyzed on a near real tims 


0aSis. 


Bie 





mx A mk WN FR © 


ema ENE eT 


“HE-NE:SCURCE ON 9,SCALE ON 
"GA-AS:SOURCE ON 8,SCALE ON 6,&CCLLIMATED CN 4": 
Pie SELECT CODE FOR PRINTER’ ,A;aev “crint” ,A 

meee, 2c -Lmt -2,L6.i,2; fmt 3,68. 1 

mune | Oca), Ri olzl¥ [5l2),Cs( 26) ,£S|[20],ES{[ 20) 

Meteo 2 | oho | oe) Be (ZU0l Ad 2,2) ,B(2z,1j) ,C( 2,1} 

bub, LATA” ,1S,22ina Y,!1 ,R:0>FewWeVeU0-1l+L:bur "I" ,IS,2 
Bece;asp "ENTER LASER USED";wait 1500 


em@e. HE-NE=1,GA-AS=6" ,h 


Mew —i7 oHE-NEB ES; 9>Cs5mp 2 


eGR-AS ES; 38>) 


(pemcoll EVAR OS OK. 5: 


beep;dsp ES&" SOURCE DATA RECCEDED?";wait 1500 


Seme.  L=-YES,0=NO" ,V 


if V=0 and w=1;stg O0;cll ‘“TRANSFER’(Q);jmp 2 
1f V=0 and W=0;srtg O;cll ‘TRANSFER’ (Q) 
beep;ent "PLOT OF SOURCE?,1=YES ,O=NC" ,R 

meek =1-cll PLOT (512,C,0) 
beep;asp ES&" SCALE CATA RECCRDED?";wait 1500 


beeprent "1=YES,G=NC" ,U 


ir U=0 and w=l;sfg O;cll “TRANSFER (7*C);imp 2 
if U=6 and W=0;sfg O;cll ‘TRANSFER’ (6+¥) 
beepsasp "wHAT IS SCALE FACTCR?";wait 1UCO 

ent “SCALE FACTOR=?;0 GETS PLOT” ,S 

ix S=6 and U#s0;U-2e0U;cll “SLOT” (512,¢,0C) 
mies—-0 ana U=O0;cll PLOT (512,0,0) 


mee s=0;jme .-3 

S> |. 

1f w=1;$+T+Qy7mp 2 
8>T=C 

ioe T,1S ,R[*) 
pos(IS,"R:") +5 

Pa l(ls(B+2,8+4] )>ri1 
memes, LT: ".)>8 
IS(B+3,5+4]+AS 

fae T-2,1S,R{ *] 
pos(IS,"R:")>B 
val(IS[B+2,5+4])>r2 
ees(is,"I1T:".)+6 
IS(B+3,Bt+4]+BS 
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me ler f—2 ,15,R[* | 


pos(I$,"R:") +5 

Ver(lso(Br2,8+4)])2rc3 

BOStIS, LI: .)*b 

I1$(B+3,B+4) cS | 

1£ r2=rl and AS=BS ang r2=r3 ana BS=CS;jmp Z 
beep;asp "ERROR IN TIME SCALE"; stp 

me slat O,1S,R[*]}eina Y;ara RY 370eA>B7C 

mom i=] tc 512 


s AtY[I)7A;:1f£ I<51;B+Y[I]->B 


if 1>463;B+Y[1})+B 

mexc is;eprt A,s 

B/1002*D;A-512*D7A 

Meemi=—) tC 52> y[1|—p=yYii|>-C+Y{1] 7c 

lt C<A/2; I-J 

Newest spr t C,J 

mele R= 1) C9256 

Minn | iwreas il 2+1)+y¥[(K+tJ-1) mced512+1] +L[K)] 
next K;ina R;ara I>kR 

Bevo POINT SPREAD FCN’ +1S:rcr 16,1$,k/*} 
beep;ent "PSF CUTPUT=l1 AND/OR CCNI" ,z 
beep;ent "PSF PLOT=l1 AND/OR CONT" ,Y 

meee =) Cll TAPLE (10) 

Mee ec) i PLOT (32,10,1) 

Moen (1 i*))+r2z2-max(1[*].)2r22 

Pee 622,023 

peep;sent “LSF CUTPUT=1 AND/CR CCNT" ,Z 
beep;ent "LSF PLOT=1 AND/OR CONT" ,Y 


cll “LSF’°:"IMAGE LINE SPREAD FCN"*IS;rcE£ 10,1S,R[*] 


maez=1,Cli TABLE (10) 
Meey=i;cil “PLOT (32,10,1) 


wet weULATICGON CP FXFCRM OF LSF" :laft 10,1$,R/[*);2ina I 


ara kR>] 

IielieR{[ 1l)s;ftor l=z tc 256 
I{[I})+R[LJ*+R[514-1I) s:next I 
pez oO) 2R( 25/7) ;recr Oye s = Ros 


Cll “FXFORM’;"FXFORM OF IMAGE LSF"*IS;rcft 10,1S,R[*] 


beep;ent "FXFORM OF LSF OUTPUT=1 AND/CR CONT" ,Z 
beep;ent "FXFORM OF LSF PLOT=1 AND/OR CCNT",Y 
iee2Z=17cll ABLE (10) 

Mayet? cl! “PLOT (32,10, 2) 

"B=DIA OF CBSCUR/DIA OF OBJ LENS IN METERS": 
"5=0.064/7.164": 

“rl=SCALE OF DATA IN MICRORADIANS,4r1>*1.60": 
"O=DIA OF OBJ LENS IN METERS;0=0.164M": 
"W=WAVELENGTH IN METERS": 





Pence Cr DIFFRACTIGN LititT PCIMA#T SPEEAD FCN": 
iné R;1it w=l;6¢.328e-7-w;jmp z 
59 .U5e- 770 
PUOo4,7, 16478: .16470;Lsrl 
B~2*D;1-D-H; 1*r1l*O*le-6/w+Z 
mom t=) EC 256 
Be (i=-LyseyY;1—E Y>sO;gtc +3 
( “AIRY (Y)-D* “AIRY (Y*B) ) °2/H”~ Z+R[T] 
Pn otO +2 
O>kK([I] 
mae asp R{l};nextel;”"OPTICS FUNCTION"=IS$ 
mere: ,1S,;R[*] 
beep;ent "OPTICS FCN PLOT=1 AND/OR CONT" ,Y 
Mmeyv=l-clil PLOT (32,11,1) 
peer cll LS >"LSF OF OPTICS FCN"+1IS;rcE— 11,15,R[*] 
beep;sent “PLOT OF LSF OF OPTICS FCN=1l1 AND/CR CONiI",Y 
Mame ecil ~ELCT (32,11,1) 
Seems; ina 1; ara Rei; L{1)]>R[1) 
fOr L=2 to z56;I[I]*R[I] eR[514-I]:next I 
Babee kRiz5>7jsccr J1,15 ,R1*)] 
cll “FXFORM’;"FXFOCRM OF LSF OF OPTICS"*IS;rcer 11,1S$,R[*] 
beep;sent "PLOT CF FAFM OF OPTICS=1 AND/OR CONTI",Y 
Maelo PLOT (256,11,2) 
Pe vecULAT PON Cr “CUOTIENT OF TWO FOURIER TRANSFORMS” : 
tie ,fS,Ri*|]sina Il;ara RI 
ieee ll, 1S,R|* ] 
Ret iL=1 tO 256 
Mil /RIiL|]s>Ri i] snext I 
aie Cr CYGIEM' sis;rct 10,135 ,k[ 7] 
beep;ent "PLOT OF MTF OF SYS=1 AND/OR CONT" ,Y 
mee iecll “PLOT “(32,10,]) 
OmerlJerlgs-rlGerz0erZler22>r2 5772477625 
~Eor [=] tc 96 
((I-1)*L)°~ (10/3) +r168 
r18+r1l727r17 
mer 1)*i) (5/3) >r19 
r19$+r20-rz0 
fer? 1] 
moi Lj) <=O;gGtCc “CC” 
(ln(R{I])27r24)* rl9tr23-r23 
r244+r257r25 
Me@e2next ©. | 
Siaemet/,6 20,021,623 ,¢625 
moe 7d )] se Zeal 1 2)]3 rl 7rAl[Z ,1j -c207A[Z,2) 
meore | l,l) 3 rZ23+C[z,1] 
inv AwA;mat A*C>B;tlt S;aprt B 
mec LCULATION OF CNSC": 
>: B(1,1]/(-21.49*145*r20*wW (.-.33333))+r22 
bist “CN S@=" ,r22>Sstp 
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ek ee 


UW) (ad Go Gd G) G) 


Stow, iar (Z 


wait 1500 

asp "PRESS CONTINUE WHEN READY"; stp 
sp aewalting= = ;walt 1500 

bur "DATA";0+2Z 

Ger 703," DATA" 

ase thanslLerring™ 

rds("CDATA") *E;1£ E=-l; jmp 0 

cies 7 U3 

asp "setting output” 

pes (iS, 's") +A; 1>N 

for I=A+2 tc 687 by 16 

aemew=1 to Je6s1£ J+tisocs; 0sru; jmp 3 
num (IS {1I+d]) 2rd 

moe y (|S | 

N+1l>N 

mext J 

next I 

remy) =) ) +B rer T=] tc Sl2;Y{I1]=-E+Y[1];next 
MicawR-ara Y>R>rCr £1,1S,Ri*];:ret 
wecreo limos tear Iiylo,R[*]samp Zz 
tee 10,1S5;,R(* | 

ina I;ara ReI;tcor L=l te 24 

feo rese ls ijl) +c 

y(I* I+U*J)>R 


7” 


Pi 


ieee aNicebhe steep; asp ENSURE EUUIF. SET UP PROPERLY” 


2* ((l-fre(R))*I[ant(k)]+tre(&)*Lf[int(R)+lj)t7¢ 


Beele git R<24;4mp. -2 
eo. {1 | 


ieee itd R-aba [>Re if cloS-rcr 11,1S,R(*];cry ;fret fl 


5 = euoue 


PHAPORM  stads9*+Nrina IL;ir tig5;lat 11,1$,R{* 


ree 10, 1S, Ri * | 

1/2° (N-1) >T 

for M=l tc N32 (N-M) +r0 

for J=0 toe 2°(M-1).-1;cll “BI °(J,P,N-1) 

cos (P¥*T)=C;sin(P*T)*(1-2*t1lg7)+P 

for T=z*r0Q*J+l1l to 2* rOG*J+r0 

RiijerisR[I+cUjer2 

meter ssi {il+rO}+r4 

mete eth aee eR Per s+ ra *C —62*P+1{ I] 

rle-rz*C-r4*P-R[It+r0]; r3-r4*C+r2* Pe Tf[I+r0] 

Mext Isnext J;csp M; next M 

Bom l=0 to 2 N-lyell B81 (1,J ,N) 

ir I-J>0;gto "5s" 

fe b=J3GtO "INC" 

ete N)>P; L{[I+l]/y i“ N)>Z 
R{[J+L]>+R[I+L] ;1[J+1] >L[I+1] 

Prieta Z> li JFL] 
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] 


;jmp 2 





PGs e tier, y (2 w)eR{[I +l} +T[itl]l) sy(z N)eL[I+1] 
wep next: | 

Peo aed 1S,n1*) saegs ret 

dqdeg;ret 

welenero2 se l>c4 

wot c=) tC ps3 

64/27p4 s2* p2-p2 

fee ( C4) FU; Ractli=cz 

int(p4) >p4 

next ZG; ret 


"AIRY":1ir pl<G;beep;asp “errcr~argument<v0";stp 


lt pl=O;l-r4;ret r4 

tees if Dl>is; jmp 2 

26+10*pl-pl° Z/3+r6; jmp :z 

SU+E1/ 2>r6 

Pe PIi<e-O64Cleriz:j jmp 2 

ie Oro cy pier 1 Z 

mex (emt tel2) int (3+e1/4))-eclz 

MOC i 2 ecOlnhg bY s7le—z6o7r6é; U>rl3>rl4 

mee Ue bee int (M7 2) pCESg LU 

Pere = eee —27 2° {h—-J)*ré/ocl—rl3erl5;r8er13 


ieee Co tia M—J—-2=b; r15>r4 


cmf£& 10; r144+2* r8*rlglOerl4;next J 

meee Ob )ser 15 

Peete loer ids £4é/rild>re4 

1x acs(r4~-r5).~-aoe (6 4*le~o) <=0; 2* r4/plerda;ret r4 
r4erS5;next M 

beee-Gse “ACCURACY NOT OBTAINEC":sweit 1500;ret ra 
Meike slat Cci,1s ,R[*] 

meeee=d EC S2erCr J=l tc 15 

mee  OLINtG. 2) ,R[1l6(L—<l)+uj];next J 

Meee OLine.s. ak | lol | snext I;ret C>2 

ener ck S22, bs,k | * } 

tA; ples;mintR(*])] )>C;max(R[*}])>+D 

peu A-,)1(B—-A) ,Bh.USIb=-A) ,C=.1(D=-C) ,p+. 05 (D-C) 
mre: LOE 

1c ¢l=512; 64>G 

me Cl=256;32+G 

1t pl=64; 5+G 

1£ pl=32;4-G 

ele B,C, 1 

mor l=— tC G by =G 

plt I(8-A) /E+A ,C,z 

mee 1 (EBAY ETA C+ (i —-C)/ i150, 2 

plt I(B-A) /E+A,C,Z 

next I 

Ber L=0 te ¢ 

plt A,I(D-C)/F+C,2Z 

teas e8-A) 7 lo, 1 (P—C) J/E+C ,2 
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Peet. Py ht eee eee eC y.Z 

next I;pen 

MS Za pst 

xa ob 

LOwel-F te 1g by ~— 1 

> plt A-.075(B-A),I(D-C)/F+C,1 

ppd 1/7 F 

next Jf 

240: Exa OQ 

fee C3=0; "DATA POINTS" +AS; 1>L 

242: if p3=2; "LINES/MICRORADIAN"+AS$;1/(2*L) +L 
Zao2 fOr I=0 tc E by G 

Pomeeolt A+(1/E-.025) (B-A) ,C-.025(D-C) , 1 

eee ID] I*L 

Beoremext Tsit e3#0 and p3#2;"MICRCRADIANS"-AS 
Pewee tk 32=2;SL 

Meee Pit .4(6-A) +A ,~~.05(D-C)+C;1leol AS 

Oe lt —.07(B=-A) +A,..5(D-C)+C;cs1z 1.z2,1,.7,90 
2503 "“NORMALIZEE INTENSITY" 2AS; 1lb1l AS 
mepweeeslZ .5,1,1.5,6;0>1 

fee Cit I ,max(kK[*]),1 

Oeeeeeeeor [=i tc pisplt I-1,R[i}];next I 

Meee cSiz i1.2,1,.7,6 

Seeeeli psrO;gtc +2 

256: beep;sent "PLOT LABEL?" ,IS 

Beemeelc .65,.SD,1;lel-"PLOT OF" 

Peeeeelt .606,.670,1;ibi I$ 

@eeeeclt .68,.840,l;lcl ES," LASER" 

meee Cen->cig -ret O-Y 

mao 29 


(Gr G? far Go GU Gt.) 
ww ®M 1700 OM & GN 
ee ee ee os ee ee 


NM po PM DW ND PS PS do 


Q4 





CATIA OP WN eH © 


Ce) 


e@ ee 8 ee e¢ ee 


"COMPLUIED SCURCE PATTERN": 


ina R; 10002A;z*L7C;trk l;:faft 10 
ror [=] tec 256 

L* (L-1) >F 

F~2/(2*C°2)>G 

me Goi3s:jmp 2 
A*texp(.-G)*R[L];1I>+K; jmp 2 

O0+>R[I] 

next I 


"COMPUTED SOURCE"+IS$5 


met 10,1S,R[*] 
Beeprent “PLCT OF CCMPUTED SOURCE=1 AND/CR CONT” ,Y 
Meevs1;cll “PLCT (32,10,G) 
"COMVERTS Peo (OC ber: 
ell “LSF 
"LSF CF COMPUTES SOURCE"+I$S 
mene 10,15 ,R[*} 
beep;ent "PLOT OF LSF OF CCMP SOURCE=1 ANCD/OR CON" ,Y 
mae iecinis PLeET (32,10, 0) 
BeALCULATION Of FXFORM OF LSF CF COMPUTED SOURCE": 
ina I;araé RI 
I{[l]>R[1] 
mer [=2 tc 256 
LiIJ]*R[I]*+R[514-I1] ;snext I 
R(c5o]*R[257};recr Oem Rae | 
Cll ‘“FXFORM 
een CF LSE OF CALC SCURCE” +IS;rcr 16,15,R/*] 
beep;sent "PLCT OF FXFM OF CALC SRCE=1 ANCD/CR CONT" ,Y 
mv l-cll PLOT (256,10,2) 
BOerLCULATICN OF PRCBUCT OF 2 FXFORMS": 
far 10,1S,k(/*|sina I;ara ReI;trk OQ 
me 1O,15,R{* } 
Or aa COZ 56 
Peek | *L[I}+R[ I} ;-next I 
mek i0,1S8,k[*] 
"TRANSFER FCN OF Poona AGecnere tonics’: 
morte ii,1S,R{[*}J;ina I;ara RI 
L{1]*R{1]>R[ 1) 
for I=Z to 256 
R[{[I] *L{[I]>R[I] *=R[(514-1} 
next I 
mez 56)]>R( 257) 
"FXFEFM OF SOURCE*ATMOS*OPTICS">=I$ 
mem 10,1S,k[ *] 
beep;ent "PLOT OF FXFM PRODUCTS=1 ANLD/OR CONT" ,Y 
mie) -C1) PEROT (32,10,2) 


aS 





"INVERSE FXFORM GIVES TARGET LSF": 

Beg @eell  FAFCRM ;cfic 7;"INVERSE FXFORM"+I$ 
mer VO, 1S, Ry | 

peep;senct "PLOT OF INV FaArM=l AND/OR CCNT" ,Y 
me YoMecwiee PLOT (32,10,0) 

"CONVERTS LSF TO PSF BY ABEL TRANSFORM": 

cll “ABEL‘;"ABEL -XFORM PSF"+IS$ 
GeiewiU,is,R[|*}];sfg 0 

beep;ent "PLOT OF ABEL XFORM=] AND/OR CONI",Y 
Meyer cli PLOT (32;,16,0)sstr 


"CALCULATES FRACTION OF POWER INSIDE CIRCLE OF RADIUS R": 


~25* 7*R[(LJ>R[1L] 

bob l=2 tc 256 

27*I*R(I]+R[L-1l]*>R[TI] 

next I 

"POWER FRCTN INSIDE CIRCLE"?IS 

mee 10,15 ,k|*} 

cpeep;ent "PLOT OF POWER=1 AND/CR CONT" ,Y 
eae cil PLOT (32,10,C) 

stp 

Shores ieee lat 1),iS ,R[*} 2 3mp 
foe ob, ilS, Rl *} 

maaeleatba: Reletor [=l tc 24 

iter aso Lei [1]} +O 

y(I*I+U*Jd)>K 

Peat LES (R))*I1( int (Rk) ] tEerc(R)*1[int(R)+4+1] )+0-vu 
Meier if R<247qnep —2 

wrI [I] 

Toei namn, ara I*Reit £lGS;rcr 11,18,R[*] ;crg 
mer pl 


NO 


eret ol 


meeoont) 2roa,5+N;ina I-it fig5;lart 11,1$,R[*] samp z 


lat 10,1S,R[(*] 

1/2°(N-1) +T 

for M=]1 tc N; Zz (N-M)+r0 

for J=0 to 2° (M-1l)-l;cli ‘BI’ (J,F,N-1) 


mpecs ( F*7T)+C;Sin(e*T) *(1-2*t1lg7) +P 


for T=2* rO*J+1l te Z2*r0*J+r0 
R[L]erl;R[I1+r0]} >r2 

peters; L{l+rC}>+r4 
rl+r2*C+r4*P+R[I] 30 34+r4*C-—r2* Pel [1] 

ieee Oar PeR[ It+cO)] -rs-r4*C+r c* P+] [i+rd j 
Mert I-next J;asp Ms next M 

Beeel=0 tc 2 N-1l:cll ‘BI (I,J ,N) 

mie r=—J>0;Ggto "BB" 

Mee i=J -gto “INC” 

R[I+l] /¥(2°N) +P; L[I+l] sy 
R[J+1l] *R[I+1l] ; I {G+l] >L[I 
P*R[J+1] ;ZeI1[J+1] 
"INC":R[I+l] /y¥(2°N) *R[I+l];1[14+l] /y(Z°N) +I [1I+1] 


(Z°N)+Z 
+1] 


S6 





TRB snext f 
aeg;ret 

PSle 20*e 2? bp l>cd 
Or sg=t EO D3 


>; p4/2-c4; z* pz-pz 


Pier me WeaeT eo E2tl +p Z 

int (p4)+p4 

next d;ret 

wwe = Pha wR; tok 10,1S ,Ri*)]>;Ri1] +N 
1.4*R[1]-1.8*R[2]h.4*R[ 3) >R[1] 
for I=2 tc 64 

R{[I]>M 
24*N+.2*M—-.6*R[I+1]*R{[I];MeNsnext I 
Bor i=] to 64 

papi, (2* 7 (ih. ty) 2-1*1))>R(1] 

for J=I+l tec 64 

R[I]+R[J] /y¥((da.1L) “2-I1*1L)+R[TI] 
next J 

Meer (i) < 20 )20sk] 1 | 
Ril)/m+R{[L);asp I;next I 

merel=95 EOC Si2;G+Riil);next I; ret 
Pebore ridin z,1s, hk! * | 
O>A;p1l-+B;min(R[*])*C;max(R[*]) +D 
SCcl A-.1(B-A) ,B8n.GS{(B-A) ,C-~.1(D-C) ,o+.05(D-C) 
Pee erortr str £lgQOsprt A,5B,C,D 

1E pl=5i2; 646 

lf pl=256; 32+G 

LZ pl=64;5+6G 

me pl=32+>4+G 

met 5,C, 1 

for [=—E tc 0 Ly -G 

plt I(E-A) /E+A,C,Z 

mice 1(S—-A) JEtA ,C+(D.-C)/ 150, 2 

met I(B-A) /Et+A,C,2z 

next I 

for I=0 to F 

ee A,1(b-C)/E+C,z 

plt A+(B-A)/150,I1(D-C)/FtC,2Z 

plt A, I(D-C)/F+tc,2z 

mext 1;rpen 

ootZ 1 .2,l,.7 

xd | 

for I=F to 0 by -l 

plt A-.G75(B-A) ,I(D-C)/Ft+C,1 

fe I/F 

next [| 

xa 0 

if p3=2; "LINES/MICRCRADIAN"+AS;1/(2*L) *L 
for L=0 to E by G 
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Meee DLlt Ar(l/E-.625) (b-&) ,C-.025(cC~-C) ,1 
me: ib) I*L 

l4o0: next I;if p3#¢2; "MICRORADIANS"+A$ 
Mawes: 10 as=2;S>L 

148; plt .4(B-A)+A,-.05(D-C)+C;1lbl AS 
meeeee le -.07 (B=-A)+A, .S(D—-Cy+CpCcsiz 1.2,1,.7,9C 
150: "NORMALIZED INTENSITY"*A$;1bl AS 
ieeeeecsiz. .5,1l,1.5,0;0+1 

fers pt I,max(R([*]),1 

ieeneecor f=) tO pl;eit I-l,R[I];next I 
Meee sSs1Z 1.2,1,./,0 

imeroit .62£,.9D,l;lbl “PLOT Or" 

Meee ple .68,.87D,1:;1lbl IS 

Meweeeett .06,.64D,1;1bl ES," LASER" 

Meee RET" spen:cig ;ret 0+Y 


71 = = «- a 
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